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ABSTRACT 

We present a detailed analysis of the Astrophysical Research Consortium 3.5m 
telescope spectrum of QSO SDSS J0838+2955. The object shows three broad 
absorption line (BAL) systems at 22,000, 13,000, and 4900 km s -1 blueshifted 
from the systemic redshift of z=2.043. Of particular interest is the lowest velocity 
system that displays absorption from low-ionization species such as Mgll, Alii, 
Sill, Sill*, Fell and Fell*. Accurate column densities were measured for all 
transitions in this lowest velocity BAL using an inhomogeneous absorber model. 
The ratio of column densities of Sill* and Fell* with respect to their ground 
states gave an electron number density of log n e (cm -3 ) = 3.75 ± 0.22 for the 
outflow. Photoionization modeling with careful regards to chemical abundances 
and the incident spectral energy distribution predicts an ionization parameter 
of log Uh = -1.93 ± 0.21 and a hydrogen column density of log N H (cm -2 ) = 
20.80 ± 0.28. This places the outflow at 3.3^}'q kpc from the central active 
galactic nucleus (AGN). Assuming that the fraction of solid angle subtended 
by the outflow is 0.2, these values yield a kinetic luminosity of (4.5^ 3 §) x 10 45 
erg s -1 , which is (1.4l ; 6 )% the bolometric luminosity of the QSO itself. Such 
large kinetic luminosity suggests that QSO outflows are a major contributor to 
AGN feedback mechanisms. 

Subject headings: galaxies: evolution, quasars: absorption lines, quasars: indi- 
vidual (SDSS J0838+2955) 
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Introduction 



Quasar outflows have been increasingly invoked as a primary feedback mechanism to 
explain the formation and evolution of supermassive black holes, th eir host gal axies, the 
surrounding intergalactic medium (IGM), and cluster cooling flows (jElvia 120061 ) . Several 
of these feedback mechanisms require that the kinetic luminosity (Ek) of the outflowing 
material would be ~5% the bolometric luminosity of the QSO itself assuming the activ e 
galactic nucleus (AGN) is accreting near the Eddington limit (jScannapieco fc Ohl 120041 ) . 
The kinetic luminosity of absorption outflows (observed as blueshifted troughs in the quasar 
spectrum) is directly proportional to the distance R of the outflow from the central AGN, 
the total hydrogen column density Njj, and the fraction Q of the solid angle subtended by 
the wind (see Section 5). These parameters must therefore be accurately determined in order 
to infer the significance of such outflows for AGN feedback. 

2002aB 



Previous attemp ts to measure E k (jWampler et al.l 119951 ; Ide Kool et al.l 12001 



Hamann et al.l 120011) had large s tatistical uncertainties (typically more than an order of 
magnitude, see iDunn et al.l 120091 ) and suffered from unquantified systematic errors due to 
the use of inadequate absorption models. The latter lead to unreliable measurements of 
column densities in the observed troughs (iVj on ), which are crucial for determining almost 
every physical aspect of the outflows. To determine Njj we use photoionization modeling 
to find the ionization equilibrium that can reproduce the measured N ion . To find R, in 
addition to the ionization equilibrium, we also need the number density of the gas which can 
be determined by comparing the ionic column densities of excited metastabl e states to their 
ground states (jWampler et al.lll995l : Ide Kool et al.ll200ll ; iKorista et al.ll2008l). We assume Q 
= 0.2 , as is the case for broad absorption line (BAL) quasars in general (IHewett & Foltz 



20031 ). We caution that Q is not directly constrained by our observations, and therefore 
express the kinetic luminosity in terms of Qq,2 (see Section 5). 

In order to determine accurate E\. in quasar outflows, we launched a multistage research 
program. We firs t developed analy sis techniques that allow us to determine reliable trough 
column densities ( lArav et al.ll2008l and references therein). We then performed an exhaustive 
search through ~50,000 QSOs with redshift z >0.2 5 and Sloan Digital Sky Survey (SDSS) 
magnitude r <19.3 in the SDSS Data Release 6 (kdelman-McCarthv et all l2008h : ~100 
objects were found that exhibited absorption systems containing transitions from excited 
metastable states of He I*, Sill*, Fell*, Fein* and/or Nill*. Of these objects, approximately 
one third exhibited only He I* (E=159,856 cm -1 ) without transitions from other energy levels 
nor other ions in metastable states that could serve as density diagnostics. Of the remaining 
~70 QSOs, ~30 objects displayed the intrinsic absorption system containing the excited 
states at velocities >1,000 km s -1 with respect to the systemic restframe of the QSO and 
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continuum flux levels >10 16 ergs 1 cm 2 A 1 in the observed frame. 

Here we investigate one of the brightest objects from this list with wide enough absorp- 
tion troughs to allow sufficiently accurate extraction of column densities with a moderate 
resolution spectrum. QSO SDSS J083817.00+295526.5 (hereafter, SDSS J0838+2955) ex- 
hibits three systems of BALs with the lowest velocity system at -4,900 km s -1 containing 
absorption from Sill and Sill* AA1527,1533 in the SDSS spectrum. With an SDSS r magni- 
tude of 17.85 and a redshift z=2.043, this object was ideal for follow-up observations because 
the optical portion of the spectrum bracketed both high ionization species of C IV and Si IV 
as well as the low ionization species of Cn, Mgll, Alii, Aim, Sill, and Fell. The plan 
of this paper is as follows. In Section [21 we detail the observations and data reduction of 
SDSS J0838+2955. Gas column and number densities are measured in Section [3] after de- 
termining the inhomogeneity of the absorber across the continuum source. Photoionization 
modeling of the outflow with careful regards to the incident spectral energy distribution 
(SED) is presented in Section HI We discuss the results for distance and kinetic luminosity 
in Section [51 For this analysis, we adopt a standa rd ACDM cosmology with Hq=70 km s _1 



Mpc 1 , fl m =0.3, and Q\=0.7 (jSpergel et al.ll2003l ). We list all transitions by their rounded 



vacuum wavelengths in the text, but use the unrounded vacuum wavelengths during analysis. 



2. Data Acquisition and Reduction 

Spectra of SDSS J0838+2955 were obtained on 2008 February 2 and 6 with the Dual 
Imaging Spectrograph (DIS) instrument on the Astrophysical Research Consortium (ARC) 
3.5 m telescope at Apache Point Observatory (APO). A total of 4.2 hr of spectra of the 
target were retrieved on the first night of observation and 3.2 hr on the second night. The 
1.5" slit and the highest resolution gratings B1200 and R1200 were used on both nights. This 
resulted in a dispersion of 0.62 A pixel -1 and 0.58 A pixel -1 in the blue and red, respectively, 
with 1.9-2.4 pixels per FWHM resolution element as measured by the calibration lamps. On 
the first night, the wavelength centers for the blue and red gratings were 4470A and 7430A 
covering the observed frames of 3850-5095A and 6825-7990A, respectively. On the second 
night, the gratings were centered at 4470A and 5745A covering the observed frames of 3850- 
5095A and 5180-6310A. This resulted in three separate regions: blue, red short, and red long 
where the blue portion received about double the exposure time of either of the red regions. 
The spectra were bias subtracted, flat-field corrected, sky subtracted, wavelength calibrated, 
extracted, atmosphere corrected with a standard star, and combined using the apall package 
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in the Imaging and Reduction Analysis Facility ( Tody 1993 1^. Considering the variable loss 



of throughput with guiding errors and cirrus clouds, an eighth degree polynomial was used 
for each region to flux calibrate the continuum in the ARC spectrum to the SDSS spectrum. 
The resulting typical FWHM resolutions were R~3300, 4500, and 5400 with continuum 
signal-to-noise ratio (S/N) per FWHM resolution element of ~105, 70, and 80 for each of 
the three regions (blue, short red, and long red). Considering the outstanding S/N but rather 
moderate resolution, the spectrum was binned to one-half the FWHM resolution element for 
analysis except for the absorption lines from Fell which were binned to one full FWHM 
resolution element to increase the S/N for these shallow troughs. The systemic redshift of 
the QSO is 2.043±0.001 based on the Siiv A1398, Aim A1859, Fein A1895, 1914, 1926, and 
Mgll A2800 emission lines observed in the ARC and SDSS spectra. 

We compare the unbinned blue region of the spectra of SDSS J0838+2955 from both 
the SDSS and ARC data in Figure [TJ Note that in this unbinned representation the ARC 
spectrum has ~3.4 times the S/N and ~1.7 times the spectral resolution compared to the 
SDSS spectrum. The three separate BAL components at 22,000, 13,000, and 4900 km s _1 
are labeled a, b, and c, while the narrow, most likely intervening systems of C IV are labeled 
1-4 where actually component 2 is a blend of C IV intervening and Si IV associated with the 
Civ system 4. Absorption in Civ and Si IV in outflow component a decreased dramatically 
between the SDSS spectrum taken in late 2003 and the ARC spectrum taken in early 2008, 
and this phenomenon is investigated further in Section [51 Note also in the figure the high 
S/N of the Sill and Sill* troughs in component c, from which an accurate number density 
can be determined. 

We normalized the continuum by fitting cubic splines across intervals of the spectrum 
that were unabsorbed. The systematic errors from normalization are negligible except in the 
case of the shallow transitions of Fe II and Fe II* where the error is a significant fraction of the 
absorbed intensity. When calculating column densities for this ion, two additional fits to the 
continuum were used to estimate the la upper and lower bounds on the continuum levels. 
We then added in quadrature the errors from normalization to the derived uncertainties of 
the column densities. 

Prior to analysis, some of the data were altered to correct contamination and other 
reduction errors. The first instance concerned the Sill* A1309 absorption trough which 
was contaminated by the red part of an intervening C IV AA1548, 1551 resonance doublet at 



4 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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z=1.528. Considering the blue component of the intervening C IV trough was uncontaminated 
with an optical depth of r ~0.18, its optical depth profile (after dividing by two because of 
the C IV doublet ratio) was used to divide out the contamination of the red component in the 
Si II* trough. The second instance was that although the Fe II A2383 line appeared centered 
at the correct location of the other troughs from component c, the next two highest strength 
lines Fell AA2600, 2344 were shifted 0.7 FWHM resolution elements to lower velocities. 
Although this much of a shift corresponds to ~6<r error based on the wavelength solution, 
underlying instrumental systematics could have caused such a shift. The two troughs were 
shifted back to match with Fell A2383 and therefore the derived column densities for Fell 
should be treated with some caution. 



3. Column and Number Densities for Component c 



Ionic column densities were determined for all the transitions covered by the ARC 
spectrum as well as Mgll which was only covered by SDSS. Atomic dat a including t ransition 
wavelen g ths, s tatistical weigh ts, and oscillato r stren gths are taken from Tayal (j2007 ) for Si II, 
Morton! (120031 ) for Fe II, and iKurucz &: Belli (119951 ) for all other ions. The adopted atomic 
parameters are displayed in Table 1. 

In order to derive the true column densities of the ionic species present in outflow com- 
ponent c, one mus t carefully consider va riations in optical depth r across the continuum 
source. Our group (lArav et al.ll2005l 120081 and references therein) showed that the apparent 
optical depth (r = -In [I] , where I is the residual intensity) method is not a good approxima- 
tion for outflow troughs (see also below). For our analysis here, we will consider the covering 
fac tor model and inh omogeneous absorber models, and use the fitting routines developed 
by lArav et al.l (120081 ). The covering factor model assumes a constant optical depth across 
the covered fraction of the continuum source and no absorption across the remainder. The 
inhomogeneous model assumes that the optical depth varies as a power-law, t(x) = r max x a , 
where the continuum source geometric span x is normalized from zero to unity and the 
power-law a parameter indicates the extent of inhomogeneity. 

We measured ionic column densities in three different ways to bracket the errors con- 
tributed by inhomogeneities of the absorber and saturation of some absorption troughs. 
Details are given below. Note that all three methods result in rather similar measurements 
of column densities (see Table 2). For seven of the ten observed species, the maximum 
column densities derived were ,$0.05 dex greater than the lower limits determined by the 
apparent optical depth method. For the three optically thick ions of Al ill, Si IV, and C IV, 
the two different inhomogeneous absorber models were able to bracket the error due to sat- 
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uration and constrain the column densities to within a factor of 2. The robustness of these 
extracted column densities lends strength to our photoionization solutions for the outflow 
(see Section 4), which in turn are used to determine the mass flux and kinetic luminosity of 
the outflow (see Section 5). 

In order to measure deviations from apparent optical depth, one must use two different 
transitions for the same species to simultaneously fit the ionic column density and the cov- 
ering factor or inhomogeneous a parameter. For component c, the Si IV AA1394, 1403 and 
Aim AA1855, 1863 resonance doublets offer a high enough S/N to constrain both param- 
eters; however, Si IV is saturated near the velocity center and Aim is contaminated from 
atmospheric emission of O I A5577 at the center (-4900 km s _1 ) of the red component and by 
an intervening absorption system at z=1.333 of Fell A2383 in the low-velocity wing (-4600 
km s _1 ) of the blue component. Nevertheless, we extracted solutions to the covering factor 
model and power-law a parameter in the inhomogeneous model for the velocity bins without 
contamination or blending, and these are shown in Figure [2j Note that the solutions to 
the covering factor are quite discrepant between the two different species of Al ill and Si IV, 
while the inhomogeneous model appears to be consistent between the two ions over the few 
velocity bins for which both ions are non-contaminated. We constructed a model for the 
power-law a parameter by smoothing the solutions with a moving boxcar average of nine 
velocity bins (dotted lines in figure). Considering the large error in the a parameter solution 
for Al ill within the high velocity wing, we used the Si IV a profile between -5900 and -5250 
km s _1 to anchor the smoothed Aim profile. Although the smoothed a profile for Aim is 
significantly lower than its actual solution at velocities -5200 to -5000 km s -1 , the derived 
total Al ill column densities differ by only ~6%. The average of the two smoothed a profiles 
is indicated by the black template line in the figure. This average template solution is used 
for reporting minimum xt values as well as ionic column densities in the figures and text, 
except for Aim and Siiv in which case their own respective smoothed solutions (dotted 
lines in figure) have been used. The average template solution is not used in Table 2 (unless 
otherwise stated), but instead the smoothed solutions for Aim and Siiv are implemented 
to bracket the errors. 

We show all of the absorption features of component c as well as the fits to the data in 
Figure [3J We integrated all column densities over the full profile of the trough (between the 
dashed lines in the figures), up to a maximum range between -6000 and -3900 km s -1 . The 
high velocity component at -6250 km s -1 is excluded because it is kinematically separated, 
seen only in the higher ionization species of Civ and Siiv and is in any case negligible, 
representing less than 2% the total true ionic column densities. Although we use only one 
power-law a parameter solution to report the column densities in the text below, we use two 
profiles to estimate the error in the column densities contributed by the uncertainty in the 
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power-law a solution. For Si IV and Al ill we used their own smoothed solutions (dotted red 
and blue lines in Figure (2], respectively) and the average template solution (black line). For 
all other ions, we compare the Si IV and Aim smoothed solutions to determine the error in 
N ion due to the uncertainty in the power-law a parameter. The measured column densities 
assuming the apparent optical depth model, the Siiv power-law a model, the Aim power- 
law a model, and the adopted solution are presented in Table 2. The adopted ionic column 
densities are derived such that they span the values and la errors of the two power-law a 
models. The details of ionic column density extractions are given henceforth. 

The ground transition of Cn A 1334. 53 is blended with the transitions from its excited 
state of E=63.4 cm' 1 at AA1335.66, 1335.71 (panel (a)). Both transitions of the excited 
state were combined into a single transition where the sum of the oscillator strengths and 
weighted wavelength transition were used. Since the electron number density log n e (cm -3 ) 
3.7 based on the Sill*/Sill ratio (see below) is above the critical density of log n e (cm -3 ) R3 
2.0 for Cn* (E=63 cm -1 ), the level population ratio of N (ground) /N (excited) is nearly 
the Boltzmann ratio which is essentially the ratio of statistical weights (1:2) for T ~ 10 4 
K. Because the ground and excited states of Cn are blended, the column density cannot 
be measured directly, and therefore a profile for the optical depth must be modeled. The 
blended profile was deconvolved by fitting Gaussian optical depth profiles to the ground and 
excited states assuming an inhomogeneous absorber model defined by the template power- 
law a parameter profile determined above (blue, red, and purple lines in panel (a)). Four 
pairs of Gaussians were used for the ground and excited state profiles, each pair with three 
free parameters (line center, optical depth, and Doppler b parameter). Minimization of the 
X 2 statistic resulted in a reduced \ 2 °f 1-7 and a column density of (8.2 ± 0.3) x 10 14 cm -2 . 

A second model for the column density profile of Cn was assumed to be the same as 
the solution for Aim (green line in panel (a)). This is motivated since Cn and Aim have 
similar ionization potentials and therefore the ratio of column densities between C II and 
Aim should be nearly constant. We multiplied the Aim column density profile by a single 
scaling factor over all velocity bins and then convolved with the power-law a parameter to 
produce the absorption profile. Two separate absorption profiles were convolved for each the 
ground and excited state, again assuming their ratio is 1:2, to produce the total absorption 
profile for the entire trough. When fitting, the error for each velocity bin is taken to be the 
error of the data for C II added in quadrature with the error in the column density solution 
of Al ill. This additional error from the Al ill column density solution is why the fit appears 
significantly worse yet with only a slightly larger reduced x 2 value of 2.2. This scaling method 
resulted in a column density of (8.8 ± 0.7) x 10 14 cm -2 , which is only 0.03 dex higher than 
the solution found from the Gaussian deconvolution solution described previously. Both the 
systematic errors that arise from the different methods employed as well as the uncertainty 
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in the a profile are accounted when reporting the adopted column density solution in Table 
2. 

The transition of Civ AA1548, 1551 is broad, blended and deep with residual intensities 
of only ~2% (panel (b)). Considering that the central portions of the Civ profile are 
saturated, the scaling of optical depth method was done with Al ill and Si IV as it was done 
for C II. We also convolved the scaled model troughs of the blue and red components of the 
Civ doublet when fitting the entire absorption profile. This is why the red component of 
the trough extends the doublet separation of ~500 km s" 1 to the right of the dashed line in 
the figure. The two different scaling fits from Aim and Siiv resulted in consistent column 
densities of (2.07 ±0.27) x 10 16 cm' 2 and (1.67 ±0.19) x 10 16 cm' 2 with a reduced x 2 of 1.6 
and 1.1, respectively. 

The transitions of Si IV AA1394, 1403 are well resolved but nearly saturated at the center 
(panel (c)). The fit using the smoothed Si IV power-law a parameter profile resulted in a 
column density of (8.4 ± 0.5) x 10 15 cm' 2 . The scaling method from the solution for Aim 
was also used, but in this case both the blue and red troughs of the doublet are fitted 
simultaneously as they are not blended and so there is no need for convolution; the light 
blue line is the scaled fit to the blue component of Si IV and the orange line is the scaled fit 
to the red component. The resulting column density of (1.12 ±0.17) x 10 16 cm" 2 is only 0.12 
dex higher than measuring Si IV on its own merits. 

The reasonable x 2 fits when scaling the optical depth profile of one ion to a different 
ion as well as the consistency in measuring column densities using the various techniques 
suggest that the ratio of column densities, even between species with different degrees of 
ionization, is fairly constant across the velocity profile. This is further evidenced by the high 
velocity wing between -5600 and -5200 km s" 1 that appear as shallow absorption features in 
the low-ionization species of Aim, Mgll, and Alii, but as deep troughs in the high ionization 
species of Civ and Siiv. Although the Civ and Si IV absorption features are ~2000 km s" 1 
wide, most of the troughs are low-column density wings while the majority of the column 
density is located at the central ~200 km s" 1 of the profile. 

For the remaining six species, Aim, Mgll, Alii, Sill, Sill* (E=287 cm" 1 ), Fell, and 
Fell* (E=385 cm" 1 ; panels (d)-(j)), we used the profile of the power-law a parameter to 
fit the column densities for each velocity bin. Although we plot only the two transitions 
AA2344, 2383 for Fell in panel (i), we used two additional lines AA2587, 2600 in the fit. The 
measured column densities assuming apparent optical depth and the various inhomogeneous 
absorber models are presented in Table 2. Note that the apparent optical depth method 
underestimates the true column densities by 0.2-0.8 dex for the optically thick, higher ion- 
ization species of Aim, Siiv, and Civ. This failure of the apparent optical depth method 



- 9 - 



Hamann et al. 


1997; 


Arav et al. 


1999 alibi. 12001. 


2002: 


Scott et al. 


2004; 


Gabel et al. 


2005a) 



To determine the electron number density of the absorber, n e , we use the column den- 
sities of the ground and first excited levels of Sill, 3s 2 3p 2 P°/ 2 an d 3s 2 3p 2 P^/ 2 & t 287 cm -1 
respectively, and of Fell, a 6 D 9 / 2 and a 6 D 7 / 2 at 385 cm -1 respectively. To ensure consistency 
when calculating number densities, the measured column densities are integrated over the 
same velocity region for each species instead of using the total column densities reported in 
Table 2. To determine t he density dependent ratio of column densities, we use the 142-level 
spectral model of Fe II of lBautista fe Pradhan (119981) and a Si H model built with atomic data 
from the CHIANTI v. 5.2 atomic database (ILandi et al.ll2006l ). both assuming an electron 
temperature T e ^10 4 K. The theoretical ratio of column densities are plotted in Figure HI 
and the measured ratio in column densities are indicated by the vertical lines. The derived 
density is simply the range in the ratio curve that spans the measured ratio. Both line ratios 
yield consistent diagnostics and we adopt the weighted average of log n e (cm -3 ) = 3.75 ± 
0.22, where the uncertainty is entirely from measurement statistics. 

With the number density now determined, the total column density in Sill is N (total) 
= A(ground) + iV(E=287 cm -1 ). To determine that of Fell we rely on our model Fell 
atom as well as the above electron density and its uncertainty to find that iV(total) = (2.5 
± 0.3) x iV(ground). These measured column densities (as well as the photoionization model 
predictions - see below) are given in Tables 3 and 4. 



4. Photoionization Analysis 

4.1. Spectral Energy Distribution 

The photon (SED) incident on the outflow is important in determining the ionization and 
thermal structures of the absorbing plasma. Because the ions used in photoionization model- 
ing have ionization potentials for destruction in the far ultraviolet (FU V), it is important to 



accura tely model this region of the SED. We first consi der the standard iMathews fe Ferland 
f jl987L hereafter MF87) SED as well as the SED used in lHamann et al.l J2OO2J, hereafter H02). 



After normalizing these SEDs near 2500 A in the rest frame of SDSS J0838+2955 to the 
flux calibrated SDSS spectrum (see Figure EJ), it becomes apparent that the observed SED 
is softer in the UV than either of these two SEDs. The slope beyond 0.4 Ryd was found to 
be a ~ -0.87 (F u oc v a \ about 0.3 dex softer than the MF87 or H02 SED. This discrepancy 
can be due to either an inherent difference in the intrinsic quasar continuum or extinction 
due to dust grains either within the outflowing gas or somewhere within the host galaxy. In 
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the latter case, grain extinction was found to diminish the rate of hydrogen ionizing photons 
by only Alog Q = 0.18±0.09 and pro duce a color excess of E(B-V)~0.03 mag, using the 
method described in lDunn et al.l (120091 ) that relies on relatively unobscured Two Micron All 
Sky Survey (2MASS) photometry data in the near IR. Considering that we do not know 
whether the deviation from the theoretical SED is due to an inherent difference in the SED, 
dust extinction that occurs within the absorber, or extinction that occurs close to the AGN, 
we simply modified the MF87 SED to match the observed UV continuum (also shown in 
Figure [5]). In addition to matching the observed slope in the UV, the UV-optical break was 
reset to 0.4 Ryd to better coincide with the SDSS spectrum, and the IR break was reset to 
5/zm to produce a more realistic heating source. We computed the bolometric luminosities 
and rate of ionizing photons Q for each of the three SEDs, and the scatter in both quantities 
are £ 0.2 dex. We consider it plausible that the actual incident FUV SED falls somewhere 
between our three adopted models. 



4.2. Photoionization Modeling 

We now examine the formation of the various ionic species observ ed in our spectrum 



For the present work, we use the photoionization modeling code Cloudy (IFerland et al.lll998f ) 
v08.00 to compute spectral models of the outflowing absorbing gas in SDSS J0838+2955. 
We start by assuming constant total hydrogen density "clouds" using Cloudy's default solar 
abundances in a plane-parallel geometry. In building a fully self-consistent photoionization 
model, we adopt a hydrogen number density n# that is 0.05 dex lower than the electron 
number density, i.e log n# (cm -3 ) = 3.70±0.22. This assumption is motivated by the 
photoionization solutions (see below) which predict the thickness of the cloud to be 0.2 
- 0.3 dex less than the column density necessary to reach the hydrogen ionization front. 
Considering the cloud contains mostly ionized hydrogen and helium, the electron number 
density in the gas region where Sill and Fell are formed is ~1.1 times higher than the 
hydrogen number density. Moreover, small variation of njj have negligible effects on the 
results of photoionization modeling. We used Cloudy's default 16-level population for Fe II, 
because for this low gas density of log n e (cm" ~ 3 ) « 3.75 at T « 10 4 K, increasing; the number 
of population levels had a negligible effect on the column densities. 

We constructed a grid of models covering the parameter space of total hydrogen column 
density Nh and ionization parameter Ujj. Here, the ionization parameter is defined as 
Uh = ^h/cuh, where (cm -2 s _1 ) is the hydrogen ionizing photon flux. The grids 
spanned log Nh (cm" 2 ) = 20.20 - 21.50 and log U H = -2 .40 to -1.40, with 0.05 dex steps 
for both parameters. The 27x21=567 individual models for a single grid were computed for 
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solar abundances only. To determine the effects of varying the abundances of silicon and 
carbon, we scaled the column densities of these species from the solar abundance grid output 
and then compared this scaled grid model to the observed column densities. This approach 
is justifiable because small variations in the silicon and carbon abundances have negligible 
effects on the ionization structure of the cloud. Thus, additional model computations are 
unnecessary. 

When comparing the ionic column densities predicted by the models to the observed 
column densities, we used two different methods to constrain Nh and Uh ■ The first approach 
entails using pairs of ionic column densities from a single element, specifically Cll and Civ 
or Sill and Siiv, to determine a solution to Nh and Uh- The strength of this method is 
that it is relatively immune to elemental abundance ratios differing from solar. Solutions 
to the Al II and Aim column density pair are not considered because there is an apparent 
problem in modeling the column dens ities of these two ions as measured in AGN outflows 



flKorista et al.ll2008l : iDunn et all 120091 1. 



The second approach in measuring N H and Uh is to co mpute the y 2 stat istic for the 



measured ionic column densities, similar to what was done in lArav et al.l (120071 ). We define 
the reduced x 2 as 



1 v-^ /log iVi 

.model log Ni 

^observed \ 



where % represents the index of one of the ionic species measured in the spectrum, a it dex is 
the measured error in dex of the ionic column densities, and v is the number of degrees of 
freedom. When minimizing x 2 , as well as determining the solution using a pair of ions, the 
grid of models is interpolated using a cubic convolution technique so that the true solution 
may be estimated. 

The solutions to Nh and Uh are shown in Figure M where the different colors represent 
the different SEDs and the various symbols represent the solutions using different constraints. 
For all SEDs, the solutions using the C II and C IV column densities (+'s) are ~0.6 dex smaller 
in Nh compared to the solutions using the Sill and Siiv column densities (O's). To study 
the effect of altering abundances, the xl statistic is minimized with the ability to vary the 
carbon abundance (d's) or the silicon abundance (A's). When computing % 2 , only four ions 
are used: Cll, Civ, Sill, and Siiv. Considering one must constrain both Nh and Uh as 
well as the chemical abundance of the single element allowed to vary, then the number of 
degrees of freedom is v — 1. We present the predicted model parameters for these solutions 
in Tables 3 and 4. Most of the solutions that use the minimization technique result in xl 
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values ~1, suggesting the fit to Nh and Uh is rather good. 

Depending on the model, the (Si/C) abundance ratio increased .1 - 0.6 dex relativ e 
to the solar ratio, as is expected in gas of super-solar metallicities ( IBallero et al.l 120081 ). 
Adopting a reasonable confidence interval that the (Si/C) abundance ratio is 0.2 - 0.5 dex 
greater than solar, then the metallicity of component c i s Z/Z Q = 1.6 - 3.1 based on the 
interpolation of the chemical abundances in lBallero et al.l (120081 . their Table 2). This super- 
solar metalli city is also consis tent with the chemical abundances determined for the outflow 
in Mrk 279 (lArav et al.l 120071 ). Even if the metallicity we re as high as Z/Z^ft ^3, the carbon 
abundance would increase by wO.l dex relative to solar (IBallero et al.l 120081 ) and therefore 
the derived solution to Nh would only be 0.1 dex less than the solution based on the carbon 
ions. 



An important conclusion drawn from the above models is that similar estimates for Nh 
and Uh are obtained despite the wide range of options in the adopted SED and chemical 
composition of the cloud. Taking the averages and standard deviations of all twelve solutions 
results in log N H (cm" 2 ) = 20.80 ± 0.28 and log U H = -1-93 ± 0.21. However, rather than 
assuming Uh to be independent from Nh, we determine a least-squares linear fit to the 
twelve solutions and adopt the correlation function to be 



log U H = -1.93 + 0.6594 x (log AVl0 2a8 cnr 2 ) ± 0.08, (2) 

which is displayed as the black parallelogram in the figure. This treatment of Uh and Nh will 
significantly reduce the errors in the derived mass of the outflow, and hence better constrain 
the mass flux and kinetic luminosity (see below). 

Figure [7] shows a representative photoionized single-slab model using the MF87 SED, 
log Uh = -1.95 and chemical abundances where carbon is decreased by 0.3 dex below solar 
while silicon is increased by 0.2 dex. We plot the ionization fraction for the species apparent 
in absorption and indicate with squares the hydrogen column densities necessary so that the 
model ionic column densities matches the observed ionic column densities for the given Uh- 
The adopted value and error in Nh are displayed as the vertical lines. Note that all the 
displayed squares are within 0.2 dex of the estimated solution for Nh- Note also that the 
ionic fractions of Cn, Mgll, Sill, and Fell rise rapi dly approaching the hydrogen ionization 



front at log Nh (cm -2 ) « 21.05 in this model (see iKorista et al.ll2008l for discussion). The 
expected column density of the absorber is 0.2 - 0.3 dex less than this necessary column to 
reach the ionization front, and therefore the lower ionization species are fairly insensitive to 
constraining Nh in this regime. 



The only ion that appears to be considerably different than the observed is Fe II, where 
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the models predict the Fe II column density should be 10 times less than measured (see Tables 
3 and 4). Considering the photoionization solutions were determined in great part by the 
silicon column densities, both the iron abundance (Fe/H) and the (Fe/Si) abundance ratio 
would have to be a full order of magnitude greater than solar to explain the discrepancy in 
the Fell column densities. However, suc h an increase in the (Fe/Si) abundance ratio is not 



predicted by chemical evolution models (jBallero et al.l 120081 ) . 



In the representative slab model, the electron temperature and electron number density 
are fairly constant across the slab with deviations from uniformity that are £ 0.2 dex and 
£j 0.05 dex, respectively. The ionic electron temperatures of Si II and Fe II were 10,400 ± 500 
in the model, justifying our assumption of T e f»10 4 K when we derived the electron number 
density from the level populations of Si II and Fe II. 

In order to compute the distance of component c to the ionizing source, we must deter- 
mine whether the flux of ionizing photons arriving to the absorber may have suffered from 
attenuation from intervening gas. Our spectrum does show strong Si IV and C IV absorption 
from a higher velocity system (component b), which might arise interior to component c. No 
absorption from lower ionization species can be measured on their own merits, suggesting 
this system is more ionized and/or lower in hydrogen column density than component c. 
In order to determine the attenuation effect, we estimated from the spectrum a la upper 
limit to the Cll column density of 9.6 x 10 13 cm -2 in component b together with a lower 
limit to N(Civ) of 8.8 x 10 15 cm" 2 from the apparent optical depth of the Civ trough of 
component b. These limits yield a lower bound of log Uh > -1.95 and an upper bound of log 
Nh (cm" 2 ) < 20.4 (cm -2 ). The attenuation of ionizing photons through such a cloud with 
these parameters is bounded by Alog Q < 0.07 dex. Thus, the uncertainty in the chosen 
SED as well as the range of solutions in Uh and Nh for the system of interest (component c) 
dominates the error compared to the uncertainty of furthered attenuation by component b. 
We therefore ignore component b and assume it has no effect on the photoionization solution 
to component c. 



5. Results and Discussion 

With estimates of nn, N h , Uh, and Q in hand, the large-scale physic al properties of the 



outflo w can be calculated. From the definition of the ionization parameter (jOsterbrock fe Ferland 



20061 . Equation (14.7)), the distance of the outflow from the central AGN is 



I 4^1=3.3-^. (3) 
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This large distance with a velocity v=4900 km s _1 suggests that this outflow will continue to 
expand into the IGM, where its mass and kinetic energy will be dissipated as heat into the 
intracluster medium. Assuming the outflow is in the form of a thin partial shell, its mass is 



M = 4nfxm p ttR 2 N H = (1.9±J;|) x 10 8 fi .2 M 



(4) 



where fi ~ 1.43 is the mean atomic mass for solar abundances, m p is the mass of the proton, 
and Q is the portion of the solid angle subtended by the shell. The adopted parameter 
space for N H and Uh is used to determine the expected uncertainty in the mass. A full 



discus sion about what Q value is appropriate for these outflows is given in iDunn et al. 



(120091 ) . To summarize their argument: there are good reasons to assume that the Q of these 
low-ionization BAL outflows with excited met ast abl e states is similar to the canonical value 
for Civ BALs (Q ~ 0.2; see iHewett fc Folta 120031 ) . and that their s carcity is due to the 



selection effect that low-ionization BALs are rare in general. See also lHall et al.l (120031 ) for 
further discussion regarding this point. In order to accommodate for the uncertainty in Q, 
we leave the final answer as a function of Qq.2 = fi/0.2. The mass flux is then given by 



M = 8-Knm p nRN H v = 590± 41 ° fi . 2 M Q yr" 1 (5) 

where the adopted parameter space for Nh and Uh is used. If the solution to Nh and Uh 
were assumed to be uncorrelated, then the error in M due to photoionization modeling would 
be 0.30 dex instead of the actual 0.19 dex. Finally, the kinetic luminosity of the outflow is 



An imipVtRN hv 6 = (4.5±f;|) x 10 45 Q . 2 erg s 



(6) 



which is (1.4 



+1.1 



0.6 ^0.2, 



the bolometric luminosity of the QSO itself. The kinetic luminosity 



of this component is sufficiently large to affect its host environment and contribute to AGN 
feedback mechanisms. Indeed, this value of ~1.4% is more likely a lower limit considering 
that the other outfl ow components a and b could further contribute to AGN feedback (see 
Arav fc et al.l 120091 ). Moreover, in Seyfert galaxies the hot phase of the outflow seen in 
X-rays can carry 70%-99% of the kinetic luminosity of the outflow (IGabel et al.l l2005bl ; 
Arav et al.ll2007l ). We point out that a few earlier investigatio ns have found outflows with 
similar physical properties. Based on the lHamann et al.l (120011 ) analysis of QSO 3C 191, the 
kinetic luminosity of the outflow observed in this object is ~9xl0 43 fio.2 er g s _1 and its mass 
flux is ~310 ftn.s M yr" 1 for an outflow situated at ~28 kpc. In QSO FBQS 1044+3656, 



de Kool et al.l (120011 ) found the outflow to be ~1 kpc from the AGN with a kinetic luminosity 
of ~9xl0 44 Qq_2 erg s _1 and mass flux of ~150 fio.2 M yr _1 . 
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Given the bolometric luminosity and assuming that the black ho le is accreting near the 
Eddington limit, we find that the amount of mass being accreted is ( ISarpeterlll964l ): 



J^Bol 

e r c 2 



60 M yr" 1 



(7) 



assuming an accretion efficiency e r of 0.1. The mass flux of the outflow is ~10 times greater 
than the accretion rate, which is physically plausible assuming the outflow entrained a sub- 
stantial amount of material on its journey to ~3 kpc away from the AGN. 

Furthermore, SDSS J0838+2955 is an example of a QSO with outflows at different 
distance scales considering that the outflow contributing to component a is most likely much 
closer than component c. The argument is as follows. The significant variation in optical 
depth of component a is most likely due to a change in flux of the io nizing source o r transverse 
velocity of the outflow out of the line of site to the central AGN ( lBarlowlll994l ). Although 
the spectrum itself cannot be used to determine variations in luminosity (because the ARC 
spectrum was flux calibrated to the SDSS spectrum), several slit guider images were used to 
conduct differential photometry between the target and foreground stars of similar color in 
the field. The change in r magnitude (the color corresponding to the peak efficiency of the 
guider CCD) was determined to be less than 0.09 mag at the 2a level and thus the variation 
in absorption of the outflow was most likely due to transverse motion. 

An upper limit on the distance of this component a of the outflow can be estimated 
assuming that the transverse velocity cannot exceed the Keplerian velocity. The Keplerian 
velocity is determined from the mass of the black hole Mbh which in turn is estimated 
from a mass-luminosity relation. In more detail, the luminosity of the QSO is expected 
to be yL y ^ 3 x 10 4 ^ erg s _1 near the restframe wavelength of 5100 A. From Equation 
(9) in iPeterson et al.l (120041 ) and assuming L/L Eddington ~0.1, the mass of the central black 



hole is expected to be Mbh ~ 7 x 10 9 M , although the mass-luminosity relation is not 
well constrained in this high-mass regime. From this black hole mass, the radius of the 
UV continuum region a rising from the accretion disk can be estimated to be ~30 light- 
days (Equation (3.20) in iPetersonl 119971 ) assuming the accretion rate is near the Eddington 
limit. We will choose a lower limit on the line of sight covering factor of component a to 
be ~0.5 because the Civ absorption trough had an initial residual intensity of ~0.43. The 
transverse velocity of the outflow across half of the ~60 light-day wide continuum source in 
the restframe time interval of 4.3 yrs / 1+z = 1.4 yrs must be v > 30 light-days / 1.4 yrs 
~ 18,000 km s _1 in the restframe of the QSO. Finally, assuming the tangential velocity of 
the outflow does not exceed the Keplerian velocity, then the distance of outflow a from the 
central AGN must be R < GM BH /v 2 « 0.1 pc. 
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Conclusions 

1. Outflow component c of SDSS J0838+2955 has a velocity of 4900 km s _1 with respect 
to the rest frame of the QSO and exhibits absorption from ground and excited metastable 
states of Sill and Fell. From the ratios of the column densities of the excited to ground 
states, the number density of the outflow was determined to be log n e (cm -3 ) = 3.75 ± 0.22. 

2. The solutions to the Aim and Siiv absorption profiles of component c indicate 
inhomogeneity of the absorber across the continuum source. The apparent optical depth 
method underestimates the true column density by 0.2-0.8 dex for the optically thick species 
of Aim, Si iv, and Civ. 

3. The uncertainty in the incident SED and chemical abundances of the outflow resulted 
in only slight variations in the derived ionization structure of the cloud. The total hydrogen 
column density for component c was determined to be log Nh (cm -2 ) = 20.80 ± 0.28, and 
the ionization parameter was found to be log Uh = -1.93 ± 0.21. The derived column density 
is 0.2 - 0.3 dex smaller than the column density necessary to reach the hydrogen ionization 
front. 

4. There is a strong indication that outflow component c has super-solar metallici- 
ties, rang ing in Z/Z^= 1.6-3.1, consistent wi th other reliable abundance estimates of AGN 



outflows flGabel et al.ll2006l : lArav et all 120071 ). 



5. We measured the distance of this outflow to be 2.3-4.8 kpc with a kinetic luminosity 
of 0.8%-2.5% the bolometric luminosity of the QSO. This large kinetic luminosity suggests 
that QSO outflows can serve as a major contributor to AGN feedback. 
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Table 1. Atomic Data for the Transitions Used 



Ion 




^vac 


9low 


log f 


Ref* 


U II 


n nn 

0.00 


1334.53 


2 


n no 

-0.92 


1 


C II 


no A O 

63.42 


i oocr 'vac 

1335.70 


4 

4 


n no 

-0.92 


1 


U IV 


r\ a A 

0.00 


1548.19 


2 


-0.72 


1 


TT T 

L IV 


r\ a A 

0.00 


1550.77 


2 


1 AO 

-1.02 


1 


Mgll 


r\ a A 

0.00 


2796.36 


2 


-0.20 


1 


Mgll 


r\ a A 

0.00 


2803.54 


2 


n n 

-0.51 


1 


A 1 TT 

Al II 


a nn 

0.00 


1670.79 


1 


n on 

0.32 


1 


A 1 TTT 
Al III 


A nn 

0.00 


1854.72 


o 

I 


n o /i 


1 


A 1 TTT 
Al III 


n nn 

0.00 


1862.79 


2 


-0.54 


1 


C * TT 

Si II 


0.00 


1304.37 


2 


-1.04 


2 


Sin* 


287.24 


1309.28 


4 


-1.09 


2 


Sin 


0.00 


1526.71 


2 


-0.88 


2 


Sin* 


287.24 


1533.43 


4 


-0.89 


2 


Sin 


0.00 


1808.01 


2 


-2.60 


2 


Sin* 


287.24 


1816.93 c 


4 


-2.67 


2 


Si iv 


0.00 


1393.76 


2 


-0.27 


1 


Si iv 


0.00 


1402.77 


2 


-0.58 


1 


Fen 


0.00 


2344.21 


10 


-0.94 


3 


Fen 


0.00 


2382.77 


10 


-0.50 


3 


Fen* 


384.79 


2396.36 


8 


-0.54 


3 


Fen 


0.00 


2586.65 


10 


-1.16 


3 


Fen 


0.00 


2600.17 


10 


-0.62 


3 



a Energy level in cm 



b References: 1- lKurucz fc Belli fll995h . 2- lTavall (120071 ) and 3- lMortonl (120031 ) 
c Blend of two transitions 
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Table 2. Comparison of Column Densities Using the Apparent Optical Depth and 

Inhomogeneous Absorber Models 



Ion 


E a 


log iV(Ap) 6 


Si IV a log iV(rnh) c 


Aim a log iV(Inh) c 


Adopted^ log 


Cll & Cll* e 


& 63 


14.88 


14.92±0.03 


14.94±0.03 


14.93±0.04 


Civ 





15.66 


16.20±0.06 


16.54±0.17 


16.42±0.28 


Mgn 





13.85 


13.87±0.08 


13.89±0.08 


13.88±0.09 


Alii 





12.88 


12.89±0.09 


12.91±0.09 


12.90±0.10 


Aim 





14.38 


14.52±0.04^ 


14.60±0.05 


14.56±0.08 


Sin 





13.88 


13.89±0.04 


13.91±0.04 


13.90±0.05 


Sin* 


287 


14.07 


14.08±0.04 


14.10±0.04 


14.09±0.05 


Si iv 





15.47 


15.93±0.04 


16.24i0.16' 


16.14±0.24 


Fen 





12.65 


12.65±0.14 


12.67±0.14 


12.66±0.15 


Fen* 


385 


12.15 


12.15±0.29 


12.16±0.29 


12.15±0.30 



a Energy level in cm 1 . 

b Apparent column density in log units of cm -2 , statistical errors are smaller than errors 
in inhomogeneous absorber models. 

c Inhomogeneous absorber column density using the indicated power-law a solution in log 
units of cm -2 . 

d The adopted column densities are derived such they span the values and la errors of the 
measurements obtained using the two power-law a models, i.e. Columns 4 and 5. 

c Assuming N(E=0)/N(E=63 cm -1 ) is the statistical weight ratio of 1:2, see text. 
f Average template power-law a solution actually used. 
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Table 3. Photoionization Models for Outflow Component c in SDSS J0838+2955 



Parameter 


SED 


MF87 


mod MF87 


H02 


MF87 


mod MF87 


H02 




Constraint 


Carbon 


Carbon 


Carbon 


Silicon 


Silicon 


Silicon 




log(JV) (cm- 2 ) 


Solar 


Solar 


Solar 


Solar 


Solar 


Solar 


log(N H ) 




20.50 


20.54 


20.54 


21.13 


21.14 


21.28 


logoff) 




-2.12 


-2.16 


-2.02 


-1.68 


-1.73 


-1.57 


log(T e ) 




4.08 


4.13 


4.13 


4.08 


4.13 


4.11 


xl 




4.2 


14 


67 


9.6 


9.2 


37 


Cn 


14.93 ±0.04 


0.00 


0.00 


0.00 


0.13 


0.14 


0.32 


Civ 


16.42 ±0.28 


0.00 


0.00 


0.00 


0.78 


0.72 


0.84 


Mgn 


13.88 ±0.09 


0.23 


-0.01 


-0.20 


0.77 


0.69 


0.73 


Alii 


12.90 ±0.10 


0.24 


0.06 


-0.10 


0.39 


0.38 


0.46 


Aim 


14.56 ±0.08 


-0.79 


-1.03 


-1.13 


-0.32 


-0.43 


-0.24 


Sin 


14.33 ±0.05 


-0.11 


-0.24 


-0.56 


0.00 


0.00 


0.00 


Si iv 


16.14 ±0.24 


-0.47 


-0.50 


-0.68 


0.00 


0.00 


0.00 


Fen 


13.03 ±0.17 


-1.14 


-1.22 


-1.34 


-1.14 


-1.10 


-1.02 



Notes, rp^g g rs ^ commn ; s the measured ionic column densities while successive columns are the predicted minus measured 
ionic column densities based on the solution for that particular model. All quantities are in log representation. Solar abundances 
have been used for these six models, and the chosen SED is indicated in the top row for each model. The first three predicted 
column densities are the solutions to Co and Civ and the last three columns arc the solutions to Sill and Si IV. The listed 
temperature is the average T e across the entire slab. The xl statistic is also reported although it was never minimized to 
determine the solution to the pair of ions. The xl value is determined only from the ClI, Civ, Sin, and Si IV values. 



Table 4. Photoionization Models for Outflow Component c in SDSS J0838+2955 



Parameter 


SED 


MF87 


mod MF87 


H02 


MF87 


mod MF87 


H02 




Constraint 


C & Si 


C & Si 


C & Si 


C & Si 


C & Si 


C & Si 




log(JV) (cm- 2 ) 


C: -0.09 


C: -0.15 


C: -0.32 


Si: ±0.14 


Si: ±0.26 


Si: ±0.60 


log(AT H ) 




20.73 


20.84 


21.03 


20.63 


20.60 


20.60 


log([/„) 




-1.97 


-1.97 


-1.77 


-2.01 


-2.11 


-1.97 


log(Te) 




4.08 


4.13 


4.10 


4.09 


4.14 


4.13 


xl 




1.8 


1.9 


1.2 


1.5 


0.93 


0.23 


Cn 


14.93 ±0.04 


0.00 


-0.01 


0.00 


-0.01 


-0.01 


-0.01 


Civ 


16.42 ±0.28 


0.21 


0.22 


0.22 


0.20 


0.11 


0.08 


Mgn 


13.88 ±0.09 


0.49 


0.48 


0.58 


0.31 


0.04 


-0.18 


Alii 


12.90 ±0.10 


0.35 


0.35 


0.46 


0.22 


0.06 


-0.11 


Aim 


14.56 ±0.08 


-0.58 


-0.65 


-0.40 


-0.72 


-1.00 


-1.10 


Sin 


14.33 ±0.05 


0.00 


0.01 


0.01 


0.01 


0.01 


0.01 


Si iv 


16.14 ±0.24 


-0.27 


-0.26 


-0.18 


-0.23 


-0.20 


-0.05 


Fen 


13.03 ±0.17 


-1.07 


-1.01 


-0.93 


-1.20 


-1.25 


-1.37 



Notes. Similar to Table 3 except the solutions are determined by minimizing the xl statistic with the option to vary the 
chemical abundance of an element. The first set of three models allow the carbon abundance to decrease while the last three 
models allow the silicon abundance to increase. The amount that the chemical abundance was varied from solar is indicated in 
the third row. 
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Fig. 1. — Blue portion of the spectrum of SDSS J0838+2955 comparing the SDSS data (red) 
taken in late 2003 and the ARC data (black) taken in early 2008. The three separate outflow 
components at 22,000, 13,000, and 4900 km s _1 are labeled a, b, and c, while the narrow, 
most likely intervening systems are labeled 1-4. Note the variability in component a as well 
as the high S/N near the Sill and Sill* absorption of component c which will allow accurate 
determination of the number density for that system. 
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Fig. 2. — Solutions for the covering factor and the power-law a parameter in the inhomoge- 
neous absorber model for Aim (blue) and Si IV (red). The covering factor solution (upper 
panel) is inconsistent between the two ions while the power-law a parameter solutions (lower 
panel) appear to overlap in the regions where both ions are uncontaminated, e.g. at veloci- 
ties -4800 to -4400 km s _1 . The inhomogeneous absorber model solutions for both ions are 
smoothed (dotted line), where the high velocity wing of Si IV is used to anchor the Aim 
profile. The average of these two smoothed solutions is indicated by the template (black), 
and all three curves will be used to estimate ionic column densities and their uncertainties. 
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Fig. 3. — Data and fit to the absorption transitions observed in outflow component c. For 
each panel, the ion and wavelengths in A of transitions are given, and the reduced 
value is indicated if the x 2 statistic was minimized. The adopted inhomogeneous absorber 
model is used to fit each velocity bin for the specific species. Gaussian fittings and scaling 
methods from other ions were employed for Cn, Civ, and Siiv (panels (a)-(c)), which all 
gave consistent results. (See text for full details.) 
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Fig. 4. — Theoretical density-dependent ratio of column densities of Si II and Si II* (E =287 
cm -1 ) as well as of Fell and Fell* (E =385 cm -1 ) for T e = 10 4 K. The uncertainty in the 
density is indicated based on the measured error of the ratios of column densities. 




Fig. 5. — Three adopted SEDs normalized to the SDSS spectrum. The bolometric luminosity 
and rate of ionizing photons Q are calculated. The ionization potentials for destruction are 
indicated by vertical lines for various ions. 
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Fig. 6. — Parameter space of total hydrogen column density Nh and ionization parameter 
Uh- The three colors represent the different SEDs used: MF87 (blue), modified MF87 (red), 
and H02 (green). The different symbols are the solutions to Nh and Uh using different 
constraints and chemical abundances. The +'s are the solutions to Nh and Uh using only 
the Cll and Civ column densities assuming solar abundances. The O's are the solutions 
using only Sill and Si IV column densities assuming solar abundances. The D's are the 
solutions to the \ 2 statistic using the four ions C II, C IV, Si II, and Si IV as constraints, but 
where the abundance of carbon was allowed to decrease. The A's are the solutions to the 
X 2 statistic where the abundance of silicon was allowed to increase. 
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Fig. 7. — Representative slab model of absorber component c assuming the MF87 SED, 
log Uh = -1.95 and chemical abundances where carbon is decreased by 0.3 dex relative to 
solar abundances while silicon is increased by 0.2 dex. The squares indicate the total Nh 
necessary so that the predicted integrated ionic column densities match the measured ones. 
The vertical lines represent the adopted value (dashed) and error (dotted) of Nh for the 
given Ujf. 



